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ABSTRACT 

We have compiled a sample of 21 low redshift {z < 0.3), luminous active 
galactic nuclei (AGN) with large Balmer decrements (Ha/H/? > 7) using the 
galaxy and QSO catalogs of the Sloan Digital Sky Survey Early Data Release. 
Using this sample we attempt to determine the fraction of quasars with large 
internal absorption. We find that these AGN have Strong [O III] A5007 and 
broad Ha emission, and that starlight dominates the spectra in the blue band, 
suggesting that these objects are heavily reddened. Their narrow emission line 
ratios are similar to those of Seyfert 2 galaxies, yet the average [O III] A5007 
emission is ~10 times more luminous. Applying the empirical relation between 
the optical continuum and the Balmer line luminosity for blue quasars, we find 
that the intrinsic luminosities of these 21 objects are in the range for quasars. 
We propose that they are obscured, intermediate type quasars analogous to type 
1.8 and 1.9 Seyfert galaxies. The ratio of these optically selected, intermediate 
type quasars to type 1 quasars are found to be around 1, similar to that for 
local Seyfert galaxies. Preliminary study indicates that most of these quasars are 
hosted in early type galaxies. 

Subject headings: galaxies: active — quasars: general — galaxies: nuclei 



1. Introduction 

The unified model for Seyfert galaxies is now well-supported (see Antonucci 1993 for a 
review). The model postulates that Seyfert 1 and Seyfert 2 galaxies are physically the same, 
except that the broad line region (BLR) in Seyfert 2 galaxies is obscured from our line of 
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sight by dust in the nucleus or in the host galaxy, and that quasars are luminous versions 
of Seyfert galaxies. Obscuration by dust should also play an important role in quasars. 
However, the evidence supporting such unification for quasars is substantially weaker than 
for Seyfert galaxies. Whether or not the difference between Seyfert galaxies and QSOs is 
simply a matter of luminosity is still a hotly debated issue, for instance some researchers 
argue that there may not even exist a dusty torus in QSOs, or, if one exists, its opening 
angle is luminosity dependent (e.g., Lawrence 1991; Simpson 2003). 

On the other hand, the hard X-ray spectrum of the cosmic X-ray background (XRB) 
seems to require a large number of obscured AGN (Seyfert galaxies and quasars). During 
the past few years, deep X-ray surveys (ROSAT, Chandra and XMM-Newton) have resolved 
more than 80% of the 2-10 keV X-ray background into discrete sources (e.g., Mushotzky et al. 
2000; Hasinger et al. 2001; Rosati et al. 2002). Follow-up spectroscopic observations show 
that the sources producing the bulk of the XRB are a mixture of obscured and unobscured 
AGN, as predicted by the XRB population synthesis models (GiUi et al. 2001 and references 
therein). However, the ratio between obscured and unobscured quasars estimated according 
to the black hole mass function in the local universe is much lower than required by the XRB 
population synthesis models (Marconi & Salvati 2002; Yu & Tremaine 2002). 

The luminous analogs of Seyfert 2 galaxies predicted by the unification model are hard 
to select and identify for the following reasons: First, they are much rarer than Seyfert 2 
galaxies, thus a large area deep survey is required. Second, it is hard to obtain their intrinsic 
luminosities because of the heavy obscuration, thus it is difiicult to distinguish them from 
Seyfert 2 galaxies. Soft X-ray surveys, such as ROSAT, can not find such objects effectively 
because the strong obscuration blocks soft X-rays. Radio surveys can be affected by free-free 
opacity associated with the obscuration, and moreover, only < 10% AGN are radio-loud 
(Krolik, 1999). Hard X-rays should be relatively unaffected by extinction, but a large-scale, 
deep survey has not yet been carried out. Though every now and then there are reports 
of the discovery of type 2 QSOs^ (e.g., Norman et al. 2002; Delia Ccca et al. 2003), the 
only systematic search for such objects covering a large area is the 2MASS red AGN survey 
(Gutri et al. 2001). 

The Sloan Digital Sky Survey (SDSS: York et al. 2000) is both an imaging and spec- 
troscopic survey. Its large sky coverage (~ 10^ deg^) makes it ideal for selecting rare objects 
such as heavily obscured quasars. In this paper we present a sample of 21 intermediate 



-'^Some of these objects may not follow the strict definition of the 'type 2' AGN, since they sometimes 
reveal weak broad Ha hnes in deeper, high S/N ratio spectra (e.g.,lE0449. 4-1823, Boyle et al. 1998; AX 
J08494-I-4454, Akiyama et al. 2002). They are, however, heavily obscured AGN. 



-3- 



type quasars selected from the SDSS EDR (Early Data Release). They are selected to have 
strong, broad Ha emission lines, with a large Balmer decrement for the broad component 
(Ha/H/? > 7), suggesting strong, partial obscuration/reddening. Their intrinsic luminosities 
are estimated from their reddening-corrected broad Ha emission. Their nearness {z < 0.3, 
with a median value of z = 0.163) makes them relatively easy for follow-up observations. 
These objects are important for three reasons: First, they can be used to study the nature of 
the obscuring material in quasars. Second, the relative number of such objects may provide 
us with some clues to the abundance of type 2 quasars, provided they are obscured by the 
same mechanism. Finally, they can be used to unravel the relationship between the nuclear 
activity and the character of the host galaxy. Usually, the study of QSO host galaxies is very 
difficult, especially for luminous QSOs, because the light from the AGN itself overwhelms 
that from the host galaxy. Thus, further work on the intermediate type QSOs discussed 
in this paper, where the nuclear emission is strongly attenuated and the host galaxy has 
prominent observed features, may be able to yield less ambiguous results on the morphology, 
stellar velocity dispersion, and star formation history in these host galaxies. 

This paper is organized as follows: In Section 2 we present the details of our analysis of 
the SDSS spectra, and describe our selection criteria for intermediate type QSOs. Section 
3 presents the main properties of our sample, as well as discussions of some interesting 
individual objects. Section 4 is a short summary. Through out the paper, we assume a 
cosmology with Ho=70 km s~^ Mpc~^, Qm—0.3 and Qa—0.7. 

2. Sample construction and Data Analysis 

2.1. Construction of the parent sample 

The parent sample of our partially obscured QSOs is the low redshift AGN (QSOs and 
Seyfert galaxies) with reliable broad Ha line selected from the spectral dataset of the SDSS 
EDR. We construct the parent sample in four stages. 

First, we pick up all the objects classified by the SDSS spectroscopic pipeline as galaxies 
(SPEG_GALAXY, specGlass=2) or quasars (SPEG_QSO, specGlass=3) with z < 0.3. This 
yields 38223 objects in total (37657 as SPEG_GALAXY, 566 as SPEG_QSO). We limit 

z < 0.3 so that the ccntroid of Ha lies at < 8500 A, allowing us to reliably measure the 
Ha flux. In this flrst step we do not limit the S/N ratio, because some objects may have 
a reliable broad Ha line yet a rather weak continuum with a very low S/N ratio. For any 
object with multiple spectroscopic observations, we retain only the spectrum with highest 
S/N ratio. 
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Second, the 566 SPEC_QSO are visually examined to remove 1) objects with too many 
bad pixels or too low S/N ratio in the Ha region, and 2) objects obviously without broad 
emission lines. This stage cuts the sample to 459 SPEC_QSO. Using these we build a sub- 
sample of normal type 1 AGN (QSOs and luminous Seyfert 1.0 galaxies, hereafter called blue 
AGN) to use as a reference in the next stages. These blue AGN are selected according to the 
following criteria: 1) -D4000 ^ 1 (the definition of the 4000 A break index -D4000 introduced 
by Balogh et al. (1999) is adopted here, using bands 3850-3890 A and 4000-4100 A) , 2) 
blue color [u' — g' ^ 0.6 corresponding to the conventional criterion oiU — B < -0.4). These 
criteria produce a sub-sample of 94 objects. For these objects, we fit the continua using a 
single power-law or broken power-law. For some objects with apparent Fe II multiplets, we 
fit the continua employing the empirical optical Fe II template of Boroson & Green (1992, 
hereafter BG92) broadened by convolving a Gaussian of various widths. The continuum 
(and Fell multiplets if present) was then subtracted and the remaining emission lines were 
fitted with an appropriate line profile (the details of the emission line fitting are described in 
§ 2.3). These emission line fits yield the equivalent widths (EW) of the broad Ha component 
of the 94 objects. These range from 134.4 to 726.2 A with a mean value of 335.4. These 
EW will be used later to roughly evaluate the relative contribution of the stellar and nuclear 
components. 

In the third stage, nuclear continuum and/or starlight subtraction is performed for 
the remaining 365 SPEC_QSO and the 37657 SPEC_GALAXY. For 41 SPEC_QSO without 
prominent stellar features, power law continuum (and Fe II templates if necessary) were fitted 
in the same way as for the blue AGN. For the other SPEC_QSO, where the stellar contribution 
cannot be neglected, nuclear continuum and stellar component are fitted simultaneously (see 
§ 2.2). 

For the 37657 SPEC_GALAXY, we initially perform starlight subtraction by fitting the 
emission line free regions of the spectrum with nine synthesized templates, which will be 
discussed in detail in § 2.2. Then the emission lines of these 37657 objects are measured 
assuming Gaussian profiles. This is a preparatory fitting for subsequent reductions. More 
than one thousand candidates were selected using the criterion that the broad component of 
Ha (Full Width at Half Maximum, FWHM> 1000 km s"^) is significant (> 5a). Then we 
visually examined these candidates to exclude the objects with a false broad component of 
Ha caused by poor starlight subtraction. After these culls, the number of the candidates is 
reduced to 867. However, for some of them, the nuclear contribution is not negligible. Those 
126 objects with EW{Ha) > 30 A (~0.2 times the minimum EW(Ha) of blue AGN and 
~0.1 times the mean value) are picked up for more precise starlight and nuclear continuum 
subtraction, as described in § 2.2. 
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The main task of the fourth stage in our analysis is to fit the emission hues precisely 
using several schemes (see § 2.3) for all the selected 1232 objects (459 SPEC_QSO and 
867 SPEC_GALAXY). Based on the emission line parameters, we select out 583 objects 
(including the 94 blue AGN) as our parent sample. The criterion is: f^""''^ > 5 (^^Habroadi 
where (y^Hairoad total flux error of the broad Ha component synthesized from two parts: 

-'-) '^Habroad^ error given by the spectral line fitting procedure and 2) (JHafiroad-: the error 
introduced by template subtraction because of the mismatch between the object's spectrum 
and the model spectrum. The latter will be discussed in § 2.2. 



2.2. steirlight/continuum subtraction 

Accurate starlight subtraction for the measurement of nuclear emission lines is impor- 
tant. Usually it is implemented by fitting the non-emission-line regions of the spectrum with 
a model spectrum composed of galaxy templates which are derived from either observed star 
spectra (e.g., Kauffmann et al. 2003) or observed galaxy spectra (e.g.. Ho et al. 1997a). In 
this work we develop a "two-step" method to build the galaxy templates described in detail 
in Li et al. (2004). Here we provide a brief summary of this approach. 

First, Principal Component Analysis (PGA) is applied to an observed stellar library 
newly available (STELIB, Le Borgne 2003). The resultant star eigenspectra broadened 
to various velocity dispersions are then used to fit the high quality spectra of more than 
one thousand galaxies selected from the SDSS spectral dataset according to 6 color-color 
diagrams using 18 synthesized magnitudes, which constitute a uniform library with a full 
spectral type coverage. Based on the fitting result, we construct model spectra of the 1016 
representative galaxies with a zero velocity dispersion. Again PGA is applied to these model 
spectra and the galaxy eigenspectra with a zero velocity dispersion are obtained. After test 
with SDSS galaxies using an F-test, the first 9 galaxy eigenspectra are picked up as our 
galaxy templates. Our galaxy templates have three advantages in addition to the merits of 
the PGA technique itself: 1) they do not include any light from the AGN, since they are 
derived from the model spectra generated from pure star eigenspectra; 2) they represent 
the most prominent features of the SDSS galaxies since they are based on representative 
galaxy spectra from the SDSS; 3) they can be used to fit the stellar velocity dispersion of 
the galaxies since they themselves are not velocity-broadened. 

When used to fit the host galaxy component of the objects, these 9 templates are 
broadened by convolving with a Gaussian to match the stellar velocity dispersion of the host 
galaxy. The host galaxy component is thus modelled as the linear combination of the 9 best 
fitting broadened spectra. 
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In our sample, many objects appear to have both a host galaxy component and a nuclear 
component. In such cases, we fit the (almost) cmission-line-free regions of the spectrum with 
the 9 broadened template spectra + a reddened power-law with Eb-v varying between 
and 2.5. The latter represents the nuclear component. Here the power-law is /(A) = CA^^'^, 
taken as the unreddened nuclear continuum according to Francis (1996). Some objects have 
visible Fe II emissions which must be also removed before spectral line fitting. For those 
objects, we first subtract the starlight/continuum component without taking account of the 
Fe II emission and fit the spectral hues. Then we broaden BG92's optical Fe II template to 
the width of the broad H/3, and fit the pseudo-continuum (wavelength ranges of prominent 
emission lines other than Fell multiplets were masked out) again with the 9 broadened 
template spectra -|- the reddened power law -|- the broadened BG92's template. 

In order to estimate the mismatch of this subtraction, we fit ~ 1000 high quality spectra 
of absorption line galaxies and obtain the standard deviation of the relative error of the 
subtraction o"r- (=0.027). Using this we estimate the error of the flux of a broad line caused 
by the subtraction as the flux of the underlying continuum within the region of 4 FWHMs 
of the line multiplied by Ur- 



2.3. Emission line fitting 

We developed a systematic code to fit the emission lines. With this code, we can fit each 
spectrum using several different fitting schemes. For example, for narrow Hct, we can tie its 
width to that of another line, such as [S II] A6731 while fitting the spectra with a extremely 
broad Ha component; or we can just fix it to a designated value, such as the fitted width 
of [O HI] A5007 when [O HI] A5007 is the only reliable narrow line. The results of different 
fitting schemes are compared with each other and the one with minimum is chosen. In 
this work, the primary goal of the line fitting is to obtain accurate measurements of broad 
Ha and broad H/3, with an emphasis on the former. To achieve this goal, we must accurately 
separate the Balmer lines from the nearby narrow lines. In general, our fitting strategies are 
as follows: 

1. Initially, lines are fitted with Gaussians, one Gaussian for each narrow or broad 
component of the lines. For example, four Gaussians are used to fit the Ha -|- [N II] lines. 
The flux ratio of the [N II] doublet A6583/A6548 is fixed to the theoretical value of 2.96, 
and their profiles are assumed to be identical. Likewise the [O HI] doublet is similarly 
constrained. Usually we also require the profiles of the narrow Ha and [N II] to be the same. 

2. If the FWHM of the broad Ha line is less than 3000 km s~^, we also try a Lorentzian 
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profile to fit the broad lines. 

3. If [O III] or [S II] doublets show a complex profile (e.g., with faint, extended wings), 
double Gaussians are employed to fit each narrow line. 

4. If the broad and narrow Ha components can not be separated by the fit, the narrow 
lines in the Ha + [N II] region will be given the same profile of the narrow H/3 line if available, 
or [S II], or [O HI] in this order. 

5. Likewise, if the broad and narrow H/5 components can not be separated by the fit, 
the narrow H/5 line will be given the same profile of the narrow Ha, or [O HI], or [S II] in 
this order. 

6. If the broad H/3 line is weak or unavailable, we constrain its profile to be that of 
broad Ha line. 

Following the above steps, we have measured the lines of most objects. But exceptions 
always exist since, as we know, the profiles of emission lines may be formed in quite different 
regions. For example, in several spectra the broad Ha component shows two or more velocity 
peaks, while in some others it even has an asymmetric or irregular shape. In these cases 
we adopt multiple Gaussians to constrain the fit if the data quality is high enough, (cf., 2 
Gaussians for the broad Ha component of SDSS J101405. 89+000620. 3 as plotted in Figure 
1). In fact, all the spectra are fitted interactively and then inspected visually. The fitting of 
the Ha + [N II] region for all the 21 partially obscured QSOs are displayed in Figure 1 (the 
right panel). 



2.4. The sample of intermediate type quaseirs 

It is well established that the Balmer line luminosity is tightly correlated with the nuclear 
continuum luminosity for both quasars and Seyfert galaxies (Yee, 1980; Shuder, 1981; Ho & 
Peng, 2001). We plot in Figure 2 the Ha luminosity versus the absolute magnitude for our 
blue AGN. A linear fit shows 

logLna = -(0.340 ± 0M8)Mg + (35.17 ± 0.40) (erg s"^) (1) 

which is perfectly consistent with the result of Ho & Peng (2001) who got 

log Lhp = -(0.340 ±0.012)M^"^+ (35.11 ±0.25) (2) 

for a sample of PG quasars and Seyfert galaxies. According to equation (1), an Ha luminosity 
of log Lffai^^g s~^) ~ 42.82 mag gives Mg — —22.5 mag as the quasar selection hmit in this 
paper. 
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Figure 3 shows the distribution of broad Ha/H/? for the 94 blue AGN we selected. It 

can be seen that the scattering of the Ha/H/S ratios is rather small, with an average Hq;/H/3 
of 2.97 ± 0.36. Assuming that all the objects in the parent sample have the same intrinsic 
Balmcr decrement, H«/H/3=2.97, we calculate the internal extinction as well as the intrinsic 
luminosity of Ha for the reddened objects. The Galactic extinction curve in Fitzpatrick 
(1999) with Rv=3.1 is utilized and the intrinsic luminosity of the broad Ha component is 
estimated as: 

log U^'^ = log Lf^ + im{\og{HalH(3) - log(2.97)) (3) 

Figure 4 shows the observed luminosity of the broad Ha component versus the Ha/H/3 
ratio for all the 583 objects in the parent sample. The tilted line corresponds to our quasar 
selection limit, i.e., M^^^^^ ~ —22.5 mag or logD^^^ ~ 42.82 erg . Out of the 80 objects 
above the tilted line, we compile a partially obscured QSO sample of 21 objects with Ha/H/9 
> 7 (listed in Table 1). The Ha/H/3 > 7 criterion is set to pick up those objects with 
relatively severe obscuration, in analogy with type 1.8/1.9 Seyfert galaxies. By this criterion, 
the sample satisfies the requirement of the intrinsic nuclear luminosity 10 times greater than 
the observed one or E(B-V) > 0.75, which is reddened more severely than the dust reddened 
quasars (E(B-V) < 0.5) compiled by Richards et al. (2003) from the SDSS QSO catalog. 
The SDSS spectra of the 21 objects are shown in Figure 1 (the left panel). 
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Table 1. Parameters of the 21 partially obscured quasars 



Object 


z 


Flux 


Flux^ 


Flux 


FWHM 


y^Est 
9 


Host 






[0 III] A5007 


H/3, broad 


Ha, broad 


Ha, broad 






•^(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


003511.48-004918.1 


.1859 


206+6 


254+17 


1847+31 


4298 


-23.2 




003916.41-003232.8 


.1098 


895+9 


73+20 


1081+21 


4135 


-22.6 




004527.06+004237.7 


.1095 


2426+15 


<125 


1968+34 


3684 


<-23.6 




015521.70-004150.0 


.2699 


482+7 


145+13 


1366+39 


4547 


-21.5 


7 


020036.60+004549.3 


.1672 


278+4 


<57 


765+25 


10004 


<-23.2 


? 


023305.96+003856.4 


.2445 


140+3 


54+11 


448+20 


4555 


-22.5 




030201.22-010117.9 


.1666 


2492+15 


<91 


1864+33 


5176 


<-25.4 


? 


031142.02-005918.9 


.2815 


123+4 


47+12 


481+30 


5093 


-23.5 


E<i 


101405.89+000620.3 


.1411 


522+8 


<340 


5509+75 


20854 


<-25.7 


E'^ 


103352.59+004403.4 


.1312 


2673+18 


356+21 


3558+41 


5465 


-23.7 




113021.41+005823.0 


.1324 


1701+12 


214+10 


2552+15 


6300 


-23.8 


E^ 




9Q1 ^ 


491 -\-A 


QQ_|_-| fl 


794-1-98 




9"^ 4 


Hi 


124321.77+001537.2 


.1426 


729+8 


54+12 


760+16 


2791 


-22.8 


Ed 


135717.59+002013.0 


.2834 


239+3 


<40 


419+33 


5513 


<-23.5 


? 


141003.66+001250.2 


.1359 


378+5 


65+12 


910+18 


3728 


-22.9 




143727.81-002343.5 


.1379 


366+6 


158+11 


1639+17 


3567 


-23.0 


E'' 


165641.98+632307.7 


.1610 


385+7 


<13 


385+12 


2109 


<-21.0 


? 


170211.15+605848.1 


.1649 


514+6 


84+25 


1036+31 


10003 


-23.3 


? 


170601.87+601732.4 


.1303 


249+4 


69+17 


1390+27 


8076 


-21.2 


E^ 


171832.86+531304.7 


.1916 


169+4 


129+14 


967+26 


2637 


-22.5 


? 


233816.44+005029.7 


.1832 


92+4 


106+11 


939+16 


3943 


-22.7 


? 



^Adopting the profile of broad Ha. 

^(1) hhmmss.ss+ddmmss.s; (3), (4) and (5) in unit of 10"^''' erg cm^^ s^^; (6) in unit of km s~^; 
(7) estimated absolute g' magnitude based on the reddening-corrected luminosity of broad H^ line; (8) 
Morphology of the host galaxy: E for early type (E, SO, Sa), S for late type. 

^Classified by visual inspection. 
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^Classified according to likelihoods and concentration index. 
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3. 



Discussion 



3.1. Testing the assumption of large, internal extinction 



We have tentatively classified 21 objects as intermediate type quasars according to their 
high intrinsic nuclear luminosities estimated by equation (3). The most remarkable property 
of these objects is their large Ha/H/3 ratios. Our basic assumption is that their intrinsic 
H«/H/9 ratio is similar to that of the "normal" QSOs, i.e. 2.97, and that the large observed 
Ha/H/3 ratio is due to relatively large internal extinction. This assumption is consistent with 
the fact that the optical continuum of almost all these objects is dominated by the stellar 
component of the host galaxy (see Figure 1). 

To test our assumption, we match the parent sample and 2MASS catalogs and find 286 
objects having counterparts within 1" in the 2MASS Point Source Catalog (PSC), without 
any extended "contamination" and with a reliable magnitude in the Kg band. Of these 
objects, the bigger the Ha/H/3 ratio, generally the redder the u' — Ks color is. All sources 
with Hq;/H/9 ratio > 7 have u' — Ks > 4.5 (Figure 5), indicating their Hq;/H/3 ratios and 
u' — Ks colors are consistent with our assumption. The absolute Ks magnitudes are well 
correlated with the luminosities of the broad Ha components of these objects (Figure 6). To 
eliminate spurious correlations due to the redshift effect, we calculate the partial correlation 
coefficient for a redshift, z, as: 



where the x, y refer to the Ha luminosity and the absolute Ks magnitude, and r^y, r^z 
and Tyz the correlation coefficients between x and x and z, and y and z. The partial 
correlation coefficient is -0.625 with a chance probability of <C0.0001 for logLna and Mxg, 
and goes up in absolute value to -0.724 for logL|£^ and A'/|^* estimated according to our 
assumption. Again, this favors our assumption of large internal extinction. 

The extinction assumption is also substantiated by the relatively high [O 111] luminosity. 
In addition to Ha luminosity, there is observational evidence that [O 111] luminosity can be 
used as a rough estimate of the intrinsic power of an AGN, even in type 2 objects (Vcron- 
Cetty & Veron 2000 and references therein). Figure 7 plots the relation between [O 111] 
luminosity and the estimated nuclear absolute magnitude for the objects with reliable [O 
HI] A5007 emission and rehable broad-hne Balmer decrements in the parent sample. The [O 
HI] luminosities are corrected for intrinsic extinction assuming an intrinsic Balmer decrement 
equal to 3.1 for the narrow components of Ha and H/3 (cf., Veilleux & Osterbrock 1987), and 
the intrinsic nuclear absolute magnitudes are estimated based on the intrinsic luminosity 



r. 



r. 



xy 



{^xz){^yz) 



(4) 
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of the broad Ha component according to equation (1) and equation (3). It is clear that 
the two are correlated, though the scatter in this relation is much larger than that for the 
correlation between the observed continuum and the Ha luminosity for the blue AGN (see 
Figure 2). The derived partial correlation coefficient is -0.578 for logL[oiii] and M*"^* with 
a chance probability of ^0.001. Hence, the correlation indicated in Figure 7 is not due to 
the redshift effect. The 87 QSOs/Seyfert Is in the BQS sample (Boroson & Green, 1992) 
as well as all the intermediate type quasars are both plotted for comparison. The [O III] 
luminosity of these intermediate quasars are comparable to that of luminous QSOs in the 
BQS sample and about 10 times more luminous than that of Seyfert galaxies, suggesting 
that they very possibly contain a quasar-like nucleus. However, from the relatively greater 
scatter of the relation between logI/[oiii] and Af*"* for both the luminous PG QSOs and our 
intermediate type QSOs, we may infer that the [O HI] luminosity is not a good calibrator 
for higher nuclear luminosity. 

The source of obscuration may be small scale dusty material near the nucleus or a large 
scale dust lane in the host galaxy. For the mean [O 111] luminosity of our intermediate 
quasars, the narrow line region (NLR) size is typically of kpc scale (Bennert et al. 2002). 
The kpc scale dust lane would also produce large reddening in the NLR, while small scale 
dusty material would not. To see which of these probabilities is most likely, we plot in 
Figure 8 the Ha/H/? ratio of the broad component versus that of the narrow component 
for 141 objects with reliably measured, narrow Ha, H/3 lines, including 9 objects which are 
intermediate type QSOs. First, it can be seen that our intermediate type quasars, together 
with most of the other partly obscured AGN, show no or moderate reddening in the NLR 
as indicated by the narrow Ha/H/9 ratio in the range 2.0 — 5.5, while some of them show 
heavy reddening in the BLR with the broad Ha/H/3 ratio > 8 and up to ~ 20. Hence, 
the absorbing material could be the dusty torus similar to that seen in Seyfert galaxies, or 
some other small scale obscuring material in the nuclear region. However a patchy, large 
scale dust lane cannot be ruled out. Second, for a rather larger number of objects, however, 
both broad and narrow lines are comparably reddened to some considerable extent (with 
the narrow and broad Ha/H/3 ratios ~ 5), indicating that the absorbing material might be 
a kpc scale dust lane or molecular disk in the host galaxy, which can play a major role in 
producing moderate, but not heavy, reddening. 



3.2. Emission lines and continua 



The broad Ha FWHMs of the 21 objects scatter widely from ~2,100 km s"^ (SDSS 
J165641.98+632307.7) to >20,000 km s"^ (SDSS J101405.89+000620.3), with a median value 
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of ~ 4500 km s^^ 

In Figure 9 we plot [O III]A5007/H/? versus [N II] /Ha and [S II] /Ha respectively for 
the 21 quasars. Most of the objects are located in the AGN region of the diagram (Veilleux 
& Osterbrock, 1987) and the distribution is similar to that of Seyfert 2 galaxies, indicating 
that these objects are similar to normal Seyfert galaxies in the properties of the NLR. Thus 

quasars, even bona fide type 2 quasars, should also bear these properties. This result is 
consistent with the recent HST observations of the NLR of bright radio-quiet PG QSOs by 
Bennert et al. (2002), who also suggested that the NLR in quasars is a scaled- up version 
of that of Seyfert galaxies. However, there is one exception. SDSS J170601. 87+601732. 4 
is located at the border between AGN and LINERs, with [O II]A3727/[0 III]A5007>1, yet 
[O I]A6300/[O III]A5007<l/3. It has an unusual spectrum and in fact is the only object in 
the parent sample that cannot be well fitted by our galaxy templates. We will discuss this 
object in detail in § 3.5. 

Our sample of intermediate type quasars can be considered to be optically selected, 
since 19 of them are optical-color selected for SDSS spectroscopic follow-up, and only 2 of 
them are selected based on ROSAT detection. We find 6 of the 21 intermediate type quasars 
have FIRST/NVSS^ counterparts within 5", and the radio detection rate is about 28±8%. 
This ratio is slightly higher than that of the 49 optically selected quasars (only 2 objects are 
selected based on ROSAT detection) at z < 0.3 in the SDSS Quasar Catalog (Schneider et 
al. 2002), which is 14%. Note an even larger radio detection rate was reported by Richards 
et al. (2001), who found that nearly half of the reddened quasars are FIRST radio sources. 
There is only one radio- loud object in each sample (Schneider et al.'s and ours). Moreover, 
each is the only source resolved by FIRST. 

3.3. The fraction of pEirtially obscured quaszirs 

In the same parent sample where we select our 21 intermediate type quasars, we also 
pick up 44 type 1 QSOs. The selection criteria we used are log LHai^ig s~^) > 42.82 (i.e., 
Mg < —22.5), and small Balmer decrements: 2.3 < Ha^^ / Hj3^^ < 4.1. The ratio of type 1 
QSOs and the intermediate type-Is is 44:21rs2.1:1. It should be pointed out that the above 
two sub-samples may not be complete, and the completeness of the samples is mainly affected 
by two factors: 1) the quality of the spectra which influenced our candidate selection and 
2) the magnitude limit that SDSS used to select targets for spectroscopic observation. The 



^Faint Images of the Radio Sky at Twenty Centimeters, Becker et al. 1995; NRAO VLA Sky Survey, 
Condon et al. 1998. 
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spectral quality does not affect tfie selection of type 1 quasars, all of which have high S/N 
spectra, however it is not negligible for the selection of intermediate type quasars. In Figure 
4 we can see that the main incompleteness comes from the measurement of the broad Hf3 
component. For many objects, due to the limited spectra quality, we can only give upper 
limits to the broad Hfi fiux. Thus our sample is conservative. A quick calculation shows that 
if we decrease the upper limits by 50%, 20 more sources would pass our selection criteria. 
We also check the possible incompleteness due to the requirement of at least a Scr detection 
of a broad Hct component. This potential selection effect is negligible, as discussed below. 
We note that over half (13 of 21) of the intermediate type QSOs are from the SDSS galaxy 
catalog with a magnitude hmit of 17.77, and that the typical S/N ratio at this magnitude 
is 17.5. We can evaluate the limiting EW of the broad Ha line with this typical S/N ratio 
and the median FWHM of the sample QSOs (4500 km s^^) corresponding to our criterion 
for reliable broad Ha line (/ > 5(t), which turn out to be 8 A. Actually, the EWs of all the 
21 objects are greater than 30 A. 

In order to correct the bias caused by the magnitude limit imposed on the spectroscopic 
target selection, for each object i we calculate -Zmax where the object would reach the magni- 
tude limit for SDSS spectroscopic follow-up. The corrected number of sources is estimated 
as E Vo.a/V^Q^, where V0.3 is the comoving volume within z = 0.3, and V^^^^ is that within 
the smaller of z^^^^ and 0.3. After the bias correction for this imposed magnitude limit, the 
number of intermediate type QSOs rises to 60.4 and the number of type 1 QSOs remains 
untouched. The bias corrected ratio of type 1 QSOs to intermediate type QSOs then be- 
comes 44:60.4Ril:1.4. To compare this with Seyfert galaxies, we pick up 12 objects from our 
intermediate type QSO sample which should be classified as type 1.8 or type 1.9 according to 
the quantitative classification introduced by Winkler (1992), and the bias corrected number 
for these type 1.8+1.9 objects is 39.1. Thus the ratio of type 1 QSOs to the type 1.8+1.9 is 
44:39.1 ~ 1:1, which is almost the same as that of Seyfert galaxies (cf., Syl.8+1.9:Syl ~ 1:1, 
Forster, 1999). But we must note that these above ratios rely on the intrinsic Ha/H/? ratio 
and the internal extinction curve we adopt for those intermediate type AGN. For example, if 
the intrinsic Hq;/H/3 ratio is set to 4.0, there are 13 intermediate type QSOs remaining and 
the bias corrected ratio of type 1 QSOs to intermediate type QSOs changes to 44:30.8^^1.4:1. 

This large fraction of obscured quasars have also been reported in many other studies. 
Through studying the red quasars in the Parkes Half-Jansky Flat-spectrum Sample, Webster 
et al. (1995) proposed that 80% of quasar population is dust-obscured and has been missed 
by optical surveys, if the dust content around radio-quiet quasars is as extensive as that 
around the radio-loud ones. Lacy et al. (2002) also suggested that the red quasar population 
could be very large through modelling the selection effects which turn out to be effective 
at removing dust- reddened quasars from magnitude- limited samples. Whiting et al. (2001), 
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however, found that ~ 40% of the Parkes sources show evidence for optical synchrotron 
emission, which is responsible for the red color of red radio-loud quasars in Webster et 
al.'s sample, as Serjeant & Rawlings (1996) argued. In that case, any missing population 
of quasars would be much smaller. But a recent study by White et al. (2003), based 
on a sample constructed by comparing the FIRST survey with a /-band survey, strongly 
supports the hypothesis that radio-loud and radio-intermediate quasars are dominated by a 
previously undetected population of red, heavily dust-obscured objects, although they do not 
know whether this kind of quasars are common in the radio-quiet population. Our sample is 
almost purely optically selected, with only one radio- loud object (see § 3.2), therefore we can 
safely conclude that dust-obscured quasars are also common in the radio-quiet population. 
As stated in §3.1, the obscuring material should be on small scales in the nucleus of the 
quasars. 

Using the SDSS quasar catalog, Richards et al. (2003) found that the redness of quasars 
whose colors cannot be fitted with a single pow-law continuum is caused by dust extinction 
instead of synchrotron emission, and moderately dust-obscured ( E(B-V) < 0.5) quasars 
account for 15% of broad-line quasars. As mentioned earlier, our sample is focused on 
heavily obscured quasars (E(B-V) > 0.75), with over half of its members from the SDSS 
galaxy catalog. Considering that Richards et al.'s sample is restricted to the SDSS quasar 
catalog, our result is marginally consistent with theirs. 

3.4. Properties of the host galeixies 

The host galaxy properties of quasars are generally difficult to determine as the stellar 
component is usually overwhelmed by the luminous quasars. High resolution ground and 
space observations are usually the only way to detect the existence of the host galaxy. Our 
sample of intermediate type quasar provides us with an excellent opportunity to study their 
host galaxies, since their nuclear emission is strongly attenuated by the obscuration. In this 
subsection, we try to draw some straightforward conclusions about the morphologies of the 
host galaxies. The study of star forming history in the host galaxies of these objects is 
deferred to another paper. 

In Figure 10, we present the SDSS images. For 5 bright objects with r' < 17 mag, we 
classify the type of host galaxy visually. For the remainder, the criteria for eUipticals are: 

• concentration index rgo/r^o > 2.5 in r* (cf. Bernardi et al. 2003); 

• the hkelihood of the de Vaucouleurs profile fit is > 0.5, and greater than the hkelihood 
of the exponential model at least by 0.1. 
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The criteria for spirals are: 

• concentration index rgo/rso < 2.5 in r*; 

• the likelihood of the exponential model fit is > 0.5, and greater than the likehhood of 
the de Vaucouleurs profile fit at least by 0.1. 

Following the above procedure, we classify 13 objects in total: 1 resides in a spiral, and 12 
reside in ellipticals (listed in Table 1). The other 8 objects are left unclassified. This is 
consistent with the result from HST observations that the host galaxies of ~2/3 QSOs are 
ellipticals (Hamilton et al. 2002). 

Some objects show possible evidence for either interacting or peculiar morphology. SDSS 
J170211. 15+605848.1 has two companions, with one at a very close projected distance (~ 
2". 5 or less than 10 kpc physical distance). SDSS J170601. 87+601732. 4 may have companion 
galaxies also. SDSS J103352. 59+004403.4 has a compact core and a very extended (50 kpc 
in diameter) low surface brightness disk viewed nearly edge on. SDSS J171832. 86+531304. 7 
has a faint spiral structure, and could be a spiral galaxy. 

3.5. Notes on individual objects 

SDSS J101405.89+000620.3: The broad Ha FWHM of this object is extremely large 
(20078+419.8 km s^^). Its broad Ha profile has a flat-top, similar to the double peaked 
line profile seen in some radio galaxies (e.g., Ostcrbrock, Koski, & Phillips 1976; Halpcrn 
& Filippenko 1988; Halpern et al. 1996). It is a 2MASS red AGN (Cutri et al. 2001) 
with a red near infrared color of J-Ks=2. 012+0. 106 mag and is bright in the near infrared 
(J=15.240±0.061, H=14.311±0.057, and K,=13.228±0.045). 

Note 19 out of the 21 intermediate type quasars are in the 2MASS Point Source Cat- 
alog, among these only 2 have J-Ks> 2, 16 objects have 1.2 < J — Kg < 2, and 1 (SDSS 
J170601.87+601732.4) has J-X, = 1.041 + 0.133. Besides SDSS J101405.89+000620.3, the 
other red object is SDSS J103352.59+004403.4, with J - K, ^ 2.073 + 0.130. This means 
that most of the partially obscured quasars would be omitted by the 2MASS red AGN survey 
{J — Kg > 2). The relatively bluer J — Kg colors may be caused by the starlight from their 
host galaxies. 

SDSS J113021.41+005823.0: This object may be classified as a FR II radio galaxy 
according to the usual classification criteria. Like SDSS J101405. 89+000620. 3, its broad Ha 
profile also has a fiat top. One interesting feature of this object is that it shows rare X- 
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shaped radio morphology with possible flows from the central radio component to the wings. 
A detailed discussion of this source was presented in Wang et al. (2003). 

SDSS J135717. 59+002013.0: This object may be taken as a heavily reddened cousin 
of the famous "post-starburst quasar" UN J1025-0040 {z = 0.634, Brotherton et al. 1999) 
located at a lower redshift {z — 0.283): its spectrum displays prominent Balmer jump and 
high-order Balmer absorption lines indicating a substantial, young stellar population. It is a 
ROSAT X-Ray source, with a fairly high luminosity in the soft X-ray band (L^ = 7.21x10'^^ 
erg s^^, uncorrected for internal extinction). It has a rather fiat ROSAT X-ray spectrum with 
photon index r=1.12]'^Q gg, indicating a heavy internal absorption. It is also fairly luminous 
in the near infrared, with = —25.9 before internal extinction correction. 

SDSS J170601. 87+601732. 4: As mentioned in § 3.2, its spectrum is the only one in the 
parent sample that cannot be well fitted even using 24 eigenspectra derived from the STELIB 
library. Its unusual spectrum is probably not be due to contamination by a foreground star, 
since the nearest star is 10" far away. However, the possibility of a foreground star sitting 
right on the galaxy can not be completely ruled out. The poor spectral fit should not be 
due to improper fiux calibration or sky subtraction either, for other spectra in the same 
plate look well-calibrated. In the light of prominent absorption lines and absorption bands 
presented in the spectrum, we propose that this spectrum is dominated by a very old and 
high metallicity stellar population, which is beyond the coverage of the STELIB library and 
thus cannot be modelled by our galaxy templates which are based on the library. As noted in 
§ 3.2, this object is LINER-like. Most LINERs are found to have predominantly old stellar 
populations (Fernandes et al. 2004; Gonzalez Delgado et al. 2004), which reinforces our 
proposal. We tend to believe that the broad Ha component in this object is real, but not 
formed by nearby absorption of neutral metals and molecules in an old stellar population, 
because: 1) its conspicuous and smooth line profile is different from that of the other smaller 
bumps blueward of it, which are likely due to nearby absorption, and 2) there is no strong 
absorption redward of the Ha emission line. Thus, SDSS J170601. 87+601732. 4 may be a 
heavily reddened, luminous version of "LINER Is" (Ho et al. 1997b). 

4. Summary and Future Prospects 

We have compiled a sample of 21 objects with relatively strong broad-Ha but very 
weak broad-H/9 emission. Assuming that the broad line is produced by photo-ionization 
process and that the large Balmer decrement is due to internal reddening, the intrinsic 
optical continuum luminosity of the nuclei, or the lower limit on it, is estimated to have 
absolute magnitude <-22'".5, implying that these are intermediate type quasars. The 
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ratio of these intermediate type quasars to type 1 quasars is similar to that of intermediate 
type Seyfert galaxies to type 1 Seyfert galaxies, indicating that quasars may be somewhat 
similar to Seyfert galaxies in structure. The narrow lines of these intermediate type quasars 
are found to be similar to those of Seyfert galaxies in line ratios, but more luminous. Our 
study has revealed that these partially obscured quasars tend to be hosted in early type 
galaxies more frequently than in late type galaxies, thus, these partially obscured quasars 
are important for the studies of the host galaxies related to morphology, stellar velocity 
dispersion and stellar population, and also the relation between the QSO phenomenon and 
the galactic evolution. A larger sample with higher S/N ratio is being constructed based on 
the observations in the SDSS Data Release 2. 

Further multi-band observations are essential to reveal the nature of these intermediate 
type quasars. Some individual objects in the sample are also very interesting and deserve 
further research. Hard X-ray observations, which can be done by XMM- Newton and/or 
Chandra (Wilkes ct al. 2002), would enable us to study the strong obscuration, and derive 
the intrinsic luminosity using an independent approach. Optical polarization spectroscopy 
would also help determine the geometry of scattering/obscuring matter inside these objects, 
using an approach similar to that carried out for infrared selected quasars (Smith et al. 
2000). 
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Fig. 1. — SDSS spectra of the intermediate type quasars and the Une fits. The vertical axis 
is flux in unit of 10~^^ erg s^^ cm~^ and the horizontal is wavelength in A. The left 
panel shows the procedure of starlight/continuum subtraction: the original spectrum (top), 
the stellar component (middle, solid line), the nuclear continuum (middle, dotted line), 
and the starhght/continuum subtracted residual (bottom). For clarity, the residual is offset 
downward by an arbitrary constant. The right panel shows the process of line parameter 
measurement. Top: original data and individual components of the fit (thin curves), final 
fit (thick curve). Bottom: the residual. 
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Fig. 1. — Continued. 
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Fig. 2. — Luminosity of the broad Ha component versus g^'-band absolute magnitude for the 
blue AGN in the parent sample. The dotted line shows a least square hnear fitting. 
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Fig. 3. — Distribution of Ha/HjS ratio of their broad components for the 94 blue AGN in 
the SDSS EDR. The dispersion of Ha/H/? ratio is rather small. 
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Fig. 4. — Luminosity of the broad Ha component versus Balmer decrement of the broad 
components of Ha and H/5 for all the AGN in the parent sample. Those which have only 
upper limit of broad H/? component are tagged with a right-pointing arrow. Blue AGN 
arc denoted solid circles. The inclined dash line corresponds to the estimated M™'^ =-22"^. 5 
after internal extinction correction. The vertical dash line denotes Ha/H/S = 7. The partially 
obscured quasars scatter in the upper-right region of the plot. 
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Fig. 5. — u'-Ks color versus the Ha;/H/3 ratio of the broad components for the 286 objects in 
the parent sample which have counterparts in the 2MASS PSC catalog without any extended 
source "contamination". Intermediate type QSOs are denoted solid squares, type 1 QSOs 
solid circles, and the remainder hollow circles. 
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Fig. 6. — Intrinsic absolute magnitude in the Ks band versus the intrinsic luminosity of the 
broad Ha component after internal extinction correction for the 286 objects, denoted as in 
Figure 5. 
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Fig. 7. — Intrinsic [O III]A5007 luminosity versus the intrinsic absolute magnitude in the g' 
band for 372 objects in the parent sample, for which [O III] flux can be accurately measured. 
The 21 intermediate type quasars (sohd square) and 87 PG QSOs (sohd triangle) are also 
plotted for comparison. Those which have only upper limit of broad H/3 component are 
tagged with a right-pointing arrow. Note that the intermediate type quasars scatter in 
almost the same area as PG QSOs. 
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Fig. 8. — Broad Ha/H/3 versus narrow Hq;/H/5 ratio of the intermediate type AGN in 
the parent sample, for which the narrow Ha and H/3 hnes can be measured accurately. 
Denotements are the same as in Figure 7. 
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Fig. 9.— The Veilleux & Osterbrock (1987) diagnostic diagram of [O III]A5007/H/3 versus [N 
II]/Hq; (the left panel) [O III]A5007/H/3 versus [N II]/Hq; (the right panel) for the 21 inter- 
mediate type quasars. The curve shows the demarcation between HII and Seyfert galaxies 
defined by Kewley et al (2001) and the horizontal straight line is [O///]A5007/if/3 = 3 
that is often used to separate Seyferts from LINERs. Note that all but 1 intermediate type 
quasars scatter in the Seyfert region. 



-38- 



fUttUMWrnu JtHESlMl-DOSmS fBVanaSHKtt^i MiS2lt94lim«S 

• # • m 

xaoms^aaansji sassDsss^sssmA JooaaiasjiiiULfi iaaiuimjiH9fi£ 



llffi«^J74DIBG03 JlflS3S2Ji^4DCM<«D4 Jtl3ini.41'IO03B£U! J]2)jnfiyUim3T'4 



JUca.iis-moi:^.] iiam^JiMDcainaji JMtixBjiisiwuicfJi ]i4}n7.«'Sd334SJ 



Fig. 10. — Composite SDSS images (from u, g, r, i and z bands) of the 21 partially obscured 
QSO candidates. The size of each image is 21"x21". 



